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INJECTION  OF  INTENSE  PCLSED-ION  BEAMS  INTO  TOKAMAKS  WITH 
APPLICATION  TO  PLASMA  HEATING  AND  CURRENT  MAINTENANCE 

1.  INTRODUCTION 


This  report  examines  the  possibility  of  injecting  intense  pulsed  ion  beams  into  a  tokamak  plasma 
to  heal  it  and  maintain  its  current.  Several  obstacles  must  be  overcome  to  accomplish  these  objectives. 
First,  the  ion  beam  must  be  produced;  second,  it  must  be  injected  into  the  toroidal  chamber;  and  third, 
it  must  deposit  its  energy  in  the  plasma.  The  crucial  problem  here  is  the  injection  [1,2].  Our  initial 
work  [ll  suggested  perpendicular  injection  into  a  steady  state  plasma.  The  present  work  concentrates 
on  a  much  easier  problem,  parallel  injection  into  the  tokamak  just  as  the  plasma  is  being  formed. 

There  are  now  several  sources  of  intense  pulsed  ion  beams,  and  the  technology  is  advancing  very 
rapidly.  For  instance  reflex  tetrodes  produce  roughly  1  MeV,  250  kA  beams  of  50  to  60  nsec  duration 
with  efficiency  of  about  55%.  Magnetically  insulated  diodes  can  produce  2  MeV,  400  kA  beams  of  85 
nsec  duration  with  80%  efficiency.  Also,  ion  beams  with  energies  in  excess  of  1  MJ  should  be  available 
very  soon  [3].  The  high  efficiencies  and  energies  make  these  beams  extremely  attractive  sources  for 
heating  plasmas. 

Several  injection  schemes  seem  to  be  viable  if  injection  takes  place  just  as  the  plasma  is  created. 
The  ion  beam  is  injected  into  a  partially  formed  plasma  and  is  trapped,  and  the  remaining  plasma  is 
formed  about  it.  This  means  either  building  up  the  density  by  fast  gas  puffing  and/or  building  up  the 
current  in  the  conventional  way  or  possibly  more  rapidly  after  the  beam  is  injected.  The  key  point  then 
is  that  the  beam  can  be  injected  only  once  and  cannot  be  used  as  an  external  power  source  in  the  con¬ 
ventional  sense. 

Whether  the  plasma  heats  up  and  then  decays,  or  else  reaches  ignition,  depends  on  the  parameters 
of  the  beam  and  plasma.  The  beam  pulse  lime  for  these  schemes  should  be  that  the  beam  (or  beams) 
just  fit  around  the  torus.  This  can  be  done  either  if  the  ion  diode  times  can  be  extended  from  about 
100  to  300  nsec,  or  if  more  than  one  injection  port  can  be  used. 

Once  the  beam  is  trapped  in  the  torus,  the  next  question  is  its  coupling  to  the  plasma.  Clearly  the 
slowing  down  time  of  the  beam  should  be  comparable  to  or  less  than  the  energy  containment  time  of 
the  plasma.  We  assume  classical  slowing  down  of  the  beam,  neoclassical  ion  losses  and  an  electron 
energy  confinement  time  which  scales  as  the  density  times  the  radius  squared  [4].  We  emphasize, 
however,  that  our  results  are  valid  only  if  this  scaling  law  still  holds  in  the  unchartered  regions  of 
parameter  space  which  we  examine. 

Our  calculations  show  that  a  tokamak  energized  with  an  intense  pulsed  ion  beam  can  achieve 
breakeven  and  the  ion  temperature  necessary  for  ignition.  For  a  small  tokamak  the  size  of  Alcator  C, 
about  1  MJ  of  beam  energy  is  needed  to  achieve  this  ion  temperature;  for  TFTR  about  40  MJ.  Our  cal¬ 
culations  show  several  specific  advantages  of  heating  with  intense  pulsed  ion  beams. 

First,  all  of  the  energy  is  in  the  plasma  at  time  r=0,  so  that  by  properly  choosing  beam  and 
plasma  parameters,  this  energy  can  be  absorbed  by  the  plasma  in  a  confinement  time  or  less.  Conven¬ 
tional  heating  schemes,  on  the  other  hand,  use  very  much  lower  power  so  that  more  injected  power  is 
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lost  from  the  plasma  as  it  slowly  heats  up.  Second,  the  beam  current  can  exceed  the  net  toroidal 
current,  allowing  for  return  current  heating.  This  heating  can  be  large  during  the  initial  heating  when 
the  plasma  temperature  is  still  low.  Thus  power  absorbed  by  the  plasma  can  be  larger  than  the  colli- 
sional  power  lost  by  the  beam.  Third,  there  is  no  limiting  plasma  density,  above  which  intense  pulsed 
ion  beams  cannot  be  used.  Fourth,  in  approaching  breakeven,  there  are  beam-plasma  reactions  as  well 
as  plasma-plasma  reactions.  Fifth,  while  most  of  the  beam  energy  is  deposited  in  electrons,  its  last  bit 
of  energy  is  deposited  directly  in  the  ions.  Since  the  electrons  are  quite  hot  at  this  time,  this  leads  to  a 
quick  boost  in  ion  temperature,  which  in  some  cases  can  lead  directly  to  ignition  temperature. 

Since  a  fixed  profile  code  is  necessary  both  to  examine  heating,  and  one  injection  scheme,  we 
begin  in  Section  II  with  the  equations  for  the  self-consistent  beam  and  plasma  response  and  show  that 
they  conserve  total  energy.  Section  111  discusses  injection  into  small-volume,  full-density  plasma  via 
Ex-B  drift  across  the  magnetic  field.  Some  preliminary  experiments  on  perpendicular  injection  have 
propagated  the  beam  across  a  field  with  almost  100‘Mi  efficiency  if  the  conditions  in  Ref.  1  are  satisfied 
[5,6].  Nearly  parallel  injection  should  be  easier  because  the  electric  fields  required  are  smaller.  Section 
IV  discusses  injection  into  a  full-volume,  low-density  plasma.  Section  V  discusses  trapping  with  a 
pulsed  stabilizing  field  or  current.  Section  VI  briefly  reviews  the  stability  of  ion  beams  in  plasmas.  Sec¬ 
tion  VII  presents  calculations  of  tokamak  heating.  Finally,  Section  Vlll  discusses  small-scale  experi¬ 
ments  that  can  be  done  to  test  some  of  the  concepts  developed  in  this  report.  In  addition  to  experi¬ 
ments  on  beam  injection  itself,  there  are  also  experiments  one  can  do  to  study  high-beta  tokamaks. 

II.  THE  EQUATIONS  FOR  BEAM  AND  PLASMA 


In  this  section  we  write  out  the  equations  for  the  self-consistent  response  of  the  beam,  plasma, 
and  poloidal  (0)  magnetic  field,  assuming  that  the  beam  has  been  injected  into  the  plasma.  The  zero 
dimensional  equations  for  the  (r  component  of  the)  beam  velocity  K,  electron  temperature  7"^ ,  and  ion 
temperature  7]  are 

dVb 

-  V^)vt,+nbeE  -  nbMt,Vi,Pbi  (1) 

=  nbM,u,(y,  -  y,)  y,  -  En^eye  -  ~  T,)  -  +  P„  -  Pr  (2) 

-  T/)  +  tihMhVbi  Uft  -  (3) 


where  n*,  and  n,  are  beam,  electron,  and  ion  number  density.  Since  and  n,  are  at  the  center  of 
the  discharge,  they  are  to  be  regarded  as  maximum,  rather  than  average  densities.  Also,  I,,  is  the  elec¬ 
tron  streaming  velocity  in  the  z  direction.  The  quantities  and  r,  are  the  electron  and  ion  energy 
confinement  times.  The  electron  confinement  time  is  taken  to  have  Alcator  scaling 

Tf  =  1.5xl0”'*nr'‘o^-  (4) 

where  r„  is  the  radius  of  the  plasma  current  channel.  The  numerical  factor  is  taken  to  agree  with  recent 
PUT  experiments.  The  ion  confinement  time  is  given  by  neoclassical  scaling 


(5) 


where 
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and 

p;=  Z,,R'^i'^Bh,r'i^  fi„(7;(ergs)/M)'/2  (8) 

where  the  units  for  T  are  specified  and  all  other  units  are  Gaussian.  The  ratio  of  ion  to  proton  mass  is 
/X.  and  \  is  the  Coulomb  logarithm.  The  collision  frequencies  are  the  beam  electron  slowing  down,  and 
electron-ion  equilibration  and  beam-ion  slowing  down  collision  frequencies,  given  by 


Vb  = 


1.7  X  10-“  n,  X 
1.6  X  10-’  n,  - 


Hi. 


1/2 


eP 


HbTP 


(eV) 


(eV), 


(9) 


whichever  is  smaller.  Here  is  the  beam  energy  in  electron  volts.  Also 


Vbi  =  9  X  10-*Z<./«<.  X 


1-h^ 


hPeP 


(eV) 


Pft?  =  3  X  10  ’^Zffrie 


(eV). 


hTP 

The  electric  field  is  gotten  from  the  z  component  of  the  electron  momentum  equation 

E  =  -IS}Ly  ( y  _  y  ) 

e  n. 


(10) 

(11) 

(12) 


where  u  is  electron-ion  momentum  exchange  collision  frequency.  The  quantities  Pa  and  Pn  are  respec¬ 
tively  the  «  particle  heating  of  electrons  (a  particle  energy  is  assumed  deposited  directly  in  the  elec¬ 
trons)  and  radiation  ioss  due  to  free  bremsstrahlung  from  the  center  of  the  plasma.  Generally  this 
latter  term  is  not  important,  but  it  can  be  important  for  high  and/or  high  density.  We  use 

P,  -  1.6  X  10-”Z^y^«/r/2(eV)  W/cm\  (13) 


The  next  thing  to  consider  is  the  equation  for  the  poloidal  magnetic  field.  This  is  given  by 


dt 


=  -  c  V  X  £|«  = 


f 


dEi 
dr  ' 


(14) 


Equation  (14)  is  inconvenient  to  work  with  because  it  has  a  strong  r  dependence.  For  simplicity,  we 
spatially  average  this  equation  in  a  way  which  conserves  total  energy.  Taking  the  dot  product  of  Eq. 
(14)  with  BgJ^,  we  find 


8  Eg 

dt  Sir 


c  cE  X 

^g,.Vx£  =  -V  — 


-  E  J. 


(15) 


The  first  term  on  ti  je  right-hand  side  of  Eq.  (15)  is  the  divergence  of  the  electromagnetic  energy  flux; 
the  second  term  denotes  an  energy  exchange  between  electromagnetic  fields  and  the  fluid.  The  quan¬ 
tity  J  is  the  total  current 


J  =  rtbeVi,  -  n^eV^. 


(16) 


Assuming  J  is  uniform  from  r  =  0  to  r  =  r,,  and  is  zero  for  <  r  <  a  where  a  is  the  liner  radius,  we 
find 


dt 


4  +  2 'm 

a 

2 

^0 

Sir 


-  Eiha^Vb  ~  n^eV^)  - 


crE^Bf 


a. 
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where  is  the  poloidal  field  at  r  =  r„.  The  {crE,BJ2)\  r  =  a  term  on  the  right-hand  side  is  the  power 
input  and  current  drive  from  the  external  circuit.  Here  we  assume  that  once  the  beam-plasma  system 
in  initialized,  there  is  no  external  power  input  or  drive.  Thus  the  equation  for  B„  is 


4-  +  2  In 

a 

^0 

B» 

-  E{nt,eVb  -  n,eV^,). 


(17) 


47r  J 

The  electron  drift  velocity  is  determined  by  Maxwell’s  equation  V  x  .  U  is 

—cB„  rik 
y  =  - * _ I — Ik 

e  '  D' 

27rr„%e 


(18) 


Equations  (1),  (2),  (3),  (12),  (17),  and  (18)  form  a  set  of  equations  for  the  time  dependence  of  the 
beam-plasma-poloidal  field  system.  The  system  conserves  energy,  since 


dt^  2 


n^Te  ^  W/T; 


+  Pa-  Pr. 


(19) 


where  is  the  total  thermal  energy  of  the  electrons  plus  ions,  plus  the  poloidal  field  energy,  plus  the 
beam  energy.  We  emphasize  once  more  that  n  this  report,  we  assume  that  after  the  beam  is  shot  in, 
there  is  no  external  power  or  current  drive  uipui  either  from  the  external  circuit  or  from  any  other 
source.  Thus,  unless  ignition  is  attained,  the  plasma  will  at  first  heat  up  and  then  decay,  that  is,  it  will 
be  transient  in  nature. 


III.  INJECTION  INTO  A  FULL-DENSITY,  LOW-VOLUME  PLASMA 

The  first  injection  scheme  we  discuss  is  the  injection  of  the  ion  beam  into  a  spatially  localized 
plasma  in  the  center  of  the  tokamak.  This  localized  plasma  is  assumed  to  carry  a  large  current,  which 
may  be  either  the  final  required  current  or  somewhat  less.  The  beam  (which  is  charge  neutralized  with 
electrons  from  the  cathode)  propagates  across  the  vacuum  region  by  setting  up  a  polarization  drift 
{1,8,91.  Then  the  E  x  B  drift  velocity  of  the  beam  will  be  nearly  equal  to  its  original  velocity  as  long  as 
(ojpft/n  »  1.  where  ojp*  and  flf*  are  respectively  the  plasma  and  cyclotron  frequency  of  the  beam 
ions.  One  great  potential  advantage  of  this  injection  scheme  is  that  a  beam-plasma  equilibrium  exists 
which  is  very  close  to  the  configuration  which  exists  at  injection. 

One  way  to  accomplish  this  injection  is  to  place  the  diode  in  a  long  guide  tube  which  joins  the 
sidewall  of  the  tokamak.  There  are  at  least  three  options,  as  shown  in  Fig.  1.  First,  the  guide  tube  may 
be  unmagnetized;  second,  the  guide  field  may  curve  away  just  before  the  entrance  to  the  tokamak;  and 
third,  the  guide  field  can  curve  and  merge  with  the  tokamak  field  like  a  bundle  diverter.  Alternately 
the  diode  may  be  placed  in  the  tokamak  (either  directly,  or  in  a  separate  chamber)  in  the  shadow  of  the 
limiter.  In  either  case,  the  beam  moves  across  the  vacuum  and  strikes  the  target  plasma.  When  it  does 
so,  it  propagates  nearly  parallel  to  the  field.  This  target  plasma  has  high  electrical  conductivity  so  that  it 
will  short  out  the  polarization  electric  field.  The  beam  will  then  be  trapped  in  the  plasma,  which  will 
serve  as  a  channel  to  guide  the  beam  around  the  torus.  Since  the  beam  does  not  carry  any  net  electrical 
current  on  injection,  the  target  plasma  must  carry  the  necessary  current  before  the  beam  is  shot  in. 
Once  the  beam  is  trapped,  the  density  is  built  up  by  gas  puffing. 

There  are  three  principal  limitations  on  this  method.  First,  the  distance  across  B  that  the  beam 
can  propagate  is  limited  by  electron  expansion  along  the  field  lines  [1]  to  roughly  d  <  (w^^/n  ,^)/), 
where  b  is  the  beam  radius.  For  our  scheme,  this  does  not  appear  to  be  a  significant  limitation. 
Second,  the  polarization  electric  field  does  not  somehow  short  circuit  (and  thereby  stop  the  beam) 
between  the  diode  and  target  plasma.  For  instance  if  the  front  end  of  the  beam  meets  the  curved  field 
while  the  rear  is  in  electrical  contact  with  the  cathode,  the  polarization  field  might  short  circuit  on  the 
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Fig.  1  —  Injection  with  an  (a)  unmagnetized,  (b)  magnetized, 
and  (c)  mostly  magnetized  guide  tube 

cathode  (9,10).  Also  the  magnetic  field  lines  in  the  tokamak  are  at  different  potentials  when  the  beam 
passes  through,  so  they  cannot  intersect  a  conductor.  Third,  the  voltage  drop  across  the  beam,  due  to 
the  polarization  field,  is  bVf,Bs\n9,  where 0  is  the  angle  between  the  field  and  beam  velocity.  Any 
injection  experiment  must  account  for  this.  Finally,  we  point  out  that  several  experiments  on  this 
injection  scheme  have  been  attempted  [5,6,11,121  and  at  least  two  [5,61  have  confirmed  the  theory  in 
Ref.  1,  and  one  [121  seems  to  work  even  better  than  predicted. 

The  experiments  in  Refs.  5  and  6  both  used  beams  with  300  and  a  configuration  like 

that  in  Fig.  1(a)  with  no  parallel  guide  field,  but  with  the  beam  perpendicular  to  the  magnetic  field.  In 
each  case  efficient  propagation  was  observed.  In  Ref.  5,  the  maximum  distance  of  propagation  was  lim¬ 
ited  to  as  predicted,  and  in  Ref.  6,  the  polarization  field  which  propagates  the  beam  was 

measured.  The  experiments  of  Ref.  11  used  a  configuration  like  that  in  Fig.  Kb)  but  only  up  to  the 
first  bend  [131.  They  found  that  a  500  kV  beam  with  =  30  propagates  across  the  bend  in  the 

field  with  more  than  50'!/ii  efficiency.  However  it  was  not  clear  whether  the  polarization  drift  or  some 
other  space  charge  effect  was  responsible  for  the  propagation 

The  University  of  California  at  Irvine  group  has  also  done  other  experiments  to  test  what  beam 
density  is  needed  to  accomplish  cross  field  injection  by  setting  up  a  polarization  field.  The\  find  that  for 
<  100,  the  electric  field  is  dominated  by  the  formation  of  a  virtual  anode  at  the  head  of  the 
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beam  (14),  whereas  for  >  300,  the  propagation  is  dominated  by  the  polarization  field  If 

beams  with  such  high  density  are  needed,  one  option  is  to  produce  them  directly  as  in  Ref.  6.  Another 
option  is  to  produce  the  beam  at  lower  density,  but  using  a  diode  with  a  rising  voltage  pulse.  The  beam 
would  then  ballistically  bunch  in  the  guide  tube  and  would  enter  the  tokamak  at  a  much  higher  density 
than  what  it  was  produced  at. 


Let  us  now  investigate  the  beam-plasma  equilibrium  which  is  likely  to  be  produced  upon  injection. 
We  assume  first  that  the  target  plasma  carries  either  the  total  tokamak  current  or  somewhat  less  (half 
the  final  current  maybe)  and  is  in  an  MUD  equilibrium.  To  be  specific,  we  assume  the  target  plasma 
has  no  pressure  so  its  equilibrium  is  force  free.  Also  we  a.ssume  that  the  relatively  low-energy  electrons 
carried  by  the  beam  immediately  mix  with  the  plasma  electrons  to  form  a  single  electron  fluid. 


The  key  to  the  new  beam-plasma  equilibrium,  and  also  the  fact  which  greatly  simplifies  it,  is  that 
current  is  frozen  into  the  fluid  on  the  beam  injection  time  scale.  To  calculate  the  new  equilibrium  one 
then  first  calculates  the  beam  orbits  in  the  existing  fields.  Secondly  one  calculates  the  response  of  the 
system  to  the  fact  that  the  current  is  now  carried  partly  by  the  beam  and  partly  by  plasma.  As  we  will 
see,  if  B  (toroidal)  >>  B  (poloidal),  the  system  hardly  responds  at  all  to  this  abrupt  change. 


We  assume  that  initially, 
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where  we  have  chosen  a  cylindrical  coordinate  system  R,(l>,z,  where  (t>  is  the  direction  of  toroidal  sym¬ 
metry.  The  flux  surfaces  are  circular,  and  the  maximum  poloidal  field  at  ro=  [(/?  -  Rq)^  +  is 
given  by  Bo  correct  to  order  rg/Ro.  The  ion  orbits  are  constrained  by 

MR  +  =  constant.  (21) 


The  singular  points  on  the  ion  orbit  are  given  by 
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If  R  Rn  is  complex,  there  are  no  singular  points,  and  therefore,  no  confined  ion  orbits.  If  /?  -  /?o  is 
real,  both  roots  are  positive,  so  that  the  orbits  are  displaced  outward  from  the  flux  surface.  The  inner 
singular  point  is  an  .v  point  and  the  outer  one  is  an  0  point  so  that  the  drift  surfaces  are  as  shown  in  Fig. 
2  The  condition  for  the  distance  between  the  two  singular  points  to  be  larger  than  ro  is 
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We  will  take  F:q.  23  as  the  condition  for  a  large  region  of  confined  ion  orbits. 


We  now  consider  the  response  of  the  plasma  to  the  presence  of  the  beam.  To  simplify  the  phy¬ 
sics.  we  replace  the  torus  with  an  equivalent  cylinder  with  coordinates  r.<t>.Z.  where  Ro(t>  =  Z.  Since 
the  plasma  is  force  free,  5/  =  iJ/  -  Ji,)B„,  where  7*  is  the  beam  current,  assumed  to  be  in  the  Z 
direction  and  J/  is  the  total  current.  Tl;us  if  B^  »  B„^  a  new  plasma  equilibrium  can  easily  form  by 
small  adjustments  in  the  poloidal  currents. 


To  conclude  we  show  that  these  adjustments  in  poloidal  currents  have  a  negligible  effect  on  the 
beam.  The  change  in  toroidal  field  produces  a  poloidal  electric  field 
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Fig  2  —  Ion  drift  orbits  in  the  poloidal  plane 


This  electric  field  gives  rise  to  a  radial  £  x  B  drift  of  the  beam,  and  the  final  radial  displacement  is 


A 


(24) 


However  since  A  Bz  «  Bz,  the  beam  displacement  is  very  small. 


To  summarize,  injecting  an  intense  pulsed  ion  beam  into  a  full-density,  low-volume  plasma  via 
the  polarization  drift  looks  very  appealing.  The  concept  itself  at  least  cleared  some  of  the  preliminary 
experimental  tests,  and  once  the  beam  is  in  the  plasma,  the  beam-plasma  system  finds  its  own  new 
equilibrium  very  easily,  without  any  help  from  an  external  circuit. 

Several  methods  could  produce  the  localized  plasma,  including  electron  cyclotron  breakdown  [15] 
and  vaporization  of  a  DT  pellet  by  a  high-power  laser.  Also  one  could  compress  either  the  toroidal  [16] 
or  vertical  field  [17]  to  draw  the  plasma  away  from  the  wall,  shoot  in  the  beam,  and  let  the  plasma  then 
expand  to  its  full  size  again.  Finally,  a  variety  of  similar  schemes  use  moving  limiters  or  suddenly 
decreased  currents  in  divertor  coils. 


IV.  BEAM  INJECTION  INTO  A  FULL-VOLUME,  LOW-DENSITY  PLASMA 

The  second  injection  scheme  we  consider  involves  beam  trapping  by  current  generation.  There 
are  two  issues  here.  First,  can  the  beam  generate  enough  current  to  stop  its  own  downward  drift,  and 
second,  can  it  generate  enough  current  to  form  a  beam-plasma  equilibrium?  The  answer  to  the  first 
question  appears  to  be  yes,  particularly  since  by  adjusting  the  vertical  field,  the  current  needed  to  stop 
the  beam  can  be  made  quite  small.  The  answer  to  the  second  question  appears  to  be  no.  Thus  once 
the  beam  is  trapped,  the  current  almost  certainly  will  have  continued  to  be  built  up  by  an  external 
source  until  the  desired  equilibrium  is  formed. 

We  consider  an  initial  full-volume,  low-density  plasma  that  carries  no  current  and  has  low  tem¬ 
perature  (in  our  calculations  we  assume  initial  temperatures  of  T,.  =  10  eV,  T,  =  5  eV').  The  beam  is 
injected  nearly  parallel  to  the  main  toroidal  field  through  an  opening  at  the  top  or  bottom  of  the  torus. 
The  density  is  low  so  that  the  return  current  velocity  is  high,  allowing  current-driven  instabilities  to  be 
excited.  Because  of  the  presence  of  the  plasma,  the  beam  is  charged  and  current  neutralized  so  the 
ions  feel  no  self  forces,  but  only  the  forces  from  the  externally  applied  toroidal  and  vertical  fields. 
Hence  each  ion  orbit  goes  principally  around  the  torus,  but  has  a  small  vertical  drift 
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where  R„  is  the  major  radius  at  injection  and  il  is  the  ion  cyclotron  freguency  in  the  toroidal  held.  We 
assume  the  drift  is  downw'ard  if  =  0. 


To  prevent  the  beam  from  just  drifting  in  the  top  and  out  the  bottom,  the  return  current  must 
decay  sufficiently  just  as  the  beam  passes  through  the  center  of  the  toroidal  chamber,  or  on  a  time  of 
order  r  =  «/>  /).  A  similar  method  of  drift  injection  has  also  been  suggested  for  electron  beams 
[18,19],  If  the  maximum  poloidal  field  is  fi„„,  one  can  show  either  by  calculating  the  drift  orbits  or  b> 
utilizing  the  constancy  of  canonical  angular  momentum  in  the  toroidal  direction,  that  the  ion  orbits  are 
confined  if  >  1  where  and  fl.,,  are  ion  cyclotron  frequencies  in  the  fields 

and  B,u  respectively.  For  Versator,  with  a  300  keV  proton  beam  with  the  drift  velocity  80'yn  canceled 
out  by  we  find  that  the  condition  for  closed  ion  orbits  is  >  300  G.  For  Alcator  C,  a  2  MeV 
tritium  beam  with  the  drift  velocity  80%  canceled  by  a  vertical  field,  requires  a  B^„  of  about  900  G. 
Thus  a  condition  for  beam  trapping  is  that  the  fields  reach  at  least  these  appropriate  values. 

The  decay  time  of  the  return  current  is  determined  in  part  by  the  background  plasma  density. 
The  drift  velocity  is  determined  by  the  vertical  field  B;„.  These  two  quantities  can  then  be  adjusted  to 
provide  optimum  conditions  for  beam  trapping. 

We  now  discuss  the  model  for  the  day  of  the  return  current.  Because  the  electron  drift  velocity  is 
so  large,  an  ion  acoustic  instability  is  excited.  This  appears  to  be  the  only  instability  which  can  act  on  a 
time  scale  sufficiently  fast  to  stop  the  beam. 


The  current-driven  ion  acoustic  instability  has  been  studied  in  many  different  laboratory  experi¬ 
ments,  and  also  in  many  different  particle  in  cell  simulations.  The  qualitative  features  of  many  of  these 
experiments  and  simulations  seem  consistent  with  ion  acoustic  turbulence  with  kXo  ~  0.5  being 
excited  to  a  fluctuating  rms  potential  of  about  0.03  <  e(i>ITf  <  0.15.  The  quasi-linear  theory  of  this 
instability  [20]  predicts  an  anomalous  electron-ion  collision  frequency  of 
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In  our  calculations  we  assume  that  Van  =  Wpp/1000  if  the  plasma  is  ion  acoustic  unstable.  Our  pro¬ 
cedure  then  is  to  add  to  r  in  Eq.  (12)  whenever  ion-acoustic  waves  are  linearly  unstable. 


The  total  plasma  heating  E  J  is  then  greatly  enhanced  by  the  anpmalous  resistance.  As  long  as  the 
ion  acoustic  waves  are  at  a  steady  state,  one  can  show  [21]  that  T^/T  ~  TJ M  >>  1.  However 
the  small  amount  of  ion  heating  plays  a  crucial  role  in  the  dynamics  because  the  condition  for  instability 
depends  on  TJT,..  If  we  neglect  the  anomalous  ion  heating,  we  find  that  much  larger  currents  can  be 
generated  and  the  trapping  can  be  greatly  improved. 


It  remains  to  determine  the  instability  threshold  for  the  ion  acoustic  waves.  This  is  actually  quite 
complicated  because  the  ion  Landau  damping  rate  depends  sensitively  on  the  ion  distribution  function. 
We  assume  that  the  ions  form  a  nonthermal  tail  on  the  distribution  function  whose  temperature  is 
roughly  the  electron  temperature.  Assuming  the  ion  heating  goes  into  producing  this  nonthermal  tail, 
the  ratio  of  the  tail  to  thermal  ion  density  is  roughly  TJT,..  If  the  ion  Landau  damping  comes  from 
this  tail,  we  find 
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(■or  this  model,  our  calculations  of  ion  temperature  are  unchanged;  the  nonthermal  tail  feature  only 
changes  the  critical  electron  drift  velocity  for  instability.  We  have  also  examined  the  effect  ol'  Maxwel¬ 
lian  ions  and  a  collision  frectuency  of  U)  -  .  The  former  choice  has  only  a  small  effect;  the  latter 

greatly  affects  the  time  scale  for  current  generation,  but  not  the  final  currents  or  temperatures. 

Calculations  were  done  for  Versator  assuming  a  beam  radius  of  5  cm.  The  unshielded  poloidal 
field  is  2400  G.  The  poloidal  field  generated  is  shown  as  a  function  of  density  in  Fig.  3.  Notice  that  at 
minimum  background  density,  nearly  half  of  the  beam  current  can  be  generated  in  the  plasma.  Also 
shown  is  the  time  at  which  the  poloidal  field  reaches  300  G  (the  field  needed  to  confine  the  beam). 
This  time  is  a  weakly  increasing  function  of  density.  For  Alcator  C  we  have  found  that  a  2  MeV,  1.2 
MA  beam  generates  1400  kG  for  a  background  density  of  3  x  lO”.  Thus  we  calculate  sufficient  current 
can  be  generated  to  at  least  stop  the  beam. 


B0  (GAUSS) 
r(n  %»c) 


Fig  .t  —  Current  generation  as  a  function  of  density 
in  Versator  for  a  5-cm  beam  radius 


We  now  turn  to  a  discussion  of  the  current  generation  experiment  on  Spac  V  done  by  Mohri  et  al. 
[19].  In  his  experiment,  500  kV.  3  cm  radius  electron  beam  was  shot  into  a  plasma  with  density 
ft  =  5  X  10'^.  The  diode  current  was  80  kA,  but  because  it  was  multiturn  injection,  the  total  current 
injected  was  about  450  kA.  He  found  about  28  kA  generated  in  about  500  nsec,  whereas  we  calculate 
about  56  kA  generated  in  500  nsec.  The  fact  that  the  time  scales  are  comparable  supports  our 
hypothesis  that  the  ion  acoustic  instability  is  responsible  for  the  anomalous  resistivity.  However  our 
calculations  show  somewhat  more  current  generated  than  what  was  measured. 

If  the  beam  generates  enough  current  to  trap  itself  but  not  enough  to  form  the  equilibrium,  the 
remaining  current  must  be  generated  externally.  One  advantage  to  using  a  low-density  background 
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plasma  is  that  the  current  can  most  likely  be  pulsed  on  rapidly  [22,23].  This  experiment  was  not  suc¬ 
cessful  as  far  as  turbulent  heating  is  concerned,  but  it  did  show  that  the  current  could  be  set  up  on  a 
microsecond  time  scale  for  a  sufficiently  low-density  background  plasma.  For  instance  with  a  =  10  cm 
and  /?  =  60  cm  and  a  background  density  of  2  x  10'^  cm  a  current  of  5  x  lO'*  A  can  be  generated  in 
1  fxsec  with  no  skin  effect.  As  the  density  increases,  the  skin  effect  becomes  more  pronounced  until  at 
5  X  10'-'  cm“^  there  is  very  little  field  penetration  [22].  The  voltages,  powers,  and  energies  involved  in 
pulsing  the  current  are  very  small  compared  to  those  involved  in  pulsing  the  beam.  Also,  it  was  shown 
that  once  the  current  is  pulsed  on,  a  normal  tokamak  discharge  can  be  produced  with  n  —  10''  cm  ' 
and  7],  ~  0.5  kV,  which  exists  in  a  steady  state  for  about  10  /xsec  and  whose  lifetime  is  limited  by  the 
external  circuit  [23].  Thus,  if  the  beam  cannot  itself  generate  enough  current  to  form  an  equilibrium, 
it  appears  that  because  the  density  is  quite  low,  the  equilibrium  can  be  produced  very  rapidly  by  an 
external  circuit. 


V.  BEAM  TRAPPING  WITH  A  PULSED  STABILIZING  OR  POLOIDAL  FIELD 

The  third  injection  scheme  involves  injection  into  a  full  volume,  charge  and  current  neutralizing 
plasma  which  carries  no  current.  The  basic  injection  process  is  the  same  as  in  Section  IV.  The  ion 
beam  comes  in  on  the  top.  Because  there  is  charge  and  current  neutrality,  the  beam  drifts  downward, 
influenced  only  by  the  externally  imposed  fields.  Assume  now  that  the  density  of  the  background 
plasma  is  high,  so  there  is  no  anomalous  resistivity.  Imagine  that  in  addition  to  the  toroidal  field,  there 
is  a  uniform  vertical  field  which  cancels  out  W’/o  of  the  downward  centrifugal  drift.  For  instance  in 
PLT,  with  a  40  kG  toroidal  field,  a  1  MeV  pri)ton  beam  and  a  1  kG  vertical  field,  the  downward  drift 
velocity  in  the  combined  toroidal  and  vertical  fields  would  be  about  3.5  x  10^  cm/ sec,  so  that  the  time 
to  drift  from  the  injector  to  the  center  of  the  toroidal  vacuum  chamber  would  be  something  over  10 
Msec.  As  the  beam  approaches  the  center,  a  stabilizing  field  would  be  pulsed  on  in  a  time  of  order  of 
this  10  Msec.  For  instance  this  stabilizing  field  might  be  a  betatron  field  [18]. 

Thus  extra  coi's  are  needed  in  the  toroidal  liner,  coils  which  can  pulse  a  field  from  say  B.  =  0.9 
B.„  to  the  betatron  field  in  perhaps  10  Msec.  This  then  will  stabilize  the  drift  near  the  center  of  the 
toroidal  chamber  and  produce  an  equilibrium  for  the  beams.  Since  this  is  a  very  small  change  in  verti¬ 
cal  field,  there  is  very  little  compression  of  the  background  plasma;  after  an  initial  slight  compression, 
the  background  plasma  will  undoubtedly  relax  back  and  contact  the  walls.  However  the  plasma  is  not  in 
equilibrium  since  its  current  is  opposite  to  the  beam  current.  Since  the  background  plasma  density  is 
high,  the  full  equilibrium  can  only  be  built  up  slowly,  probably  on  a  millisecond  time  scale.  Therefore 
this  injection  scheme  can  be  viable  only  if  the  beam-plasma  system  can  exist  without  equilibrium,  say 
in  some  complex  plasma  flow  pattern,  for  these  long  times. 

VL  REVIEW  OF  STABILITY  OF  INTENSE  BEAMS  IN  PLASMAS 


The  equilibrium  for  the  beam  plasma  system  is  simpler  than  those  worked  out  previously  [24] 
because  the  plasma  currents  <as  discussed  in  Sec.  Ill*  play  an  important  rule  in  forming  the  equilibrium 
We  now  turn  to  the  question  of  the  stability  of  the  beam  plasma  system  [25-30].  One  useful  result 
derived  by  Lovelace  [25]  for  a  large  aspect  ratio  torus  is  that  if  the  tokamak  plasma  has  no  equilibrium 
current  or  pressure,  the  .MHD  stability  of  the  beam  plasma  system  follows  from  the  MHO  stability  of 
the  plasma  itself,  if  one  makes  the  replacement  p  ^  Thus  one  has  the  Kruskal- 

Shafranov  stability  condition.  Also  one  has  a  simple  analog  to  the  Suydam  stability  condition  which  is 
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where  /3,  =  An  Br  and  q  =  iBJr)l  RBu^r).  Lovelace  also  finds  a  sufficient  stability  condition 

for  kink  modes  with  w  =  1.  0  <  n'  «  iRIa)'  to  be 
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Another  possible  MUD  instability  is  the  ballooning  mode.  If  Lovelace’s  result,  p  — 
holds  in  a  torus,  as  well  as  a  straight  cylinder,  ballooning  modes  can  be  driven  by  the  beam  energy  den¬ 
sity,  just  as  they  are  driven  by  pressure  in  ideal  MUD.  This  could  be  a  significant  problem  because  the 
beam  energies  we  deal  with  are  not  that  small  compared  to  the  magnetic  energy.  In  our  V’ersator  exam¬ 
ple,  the  local  /3,  is  about  15%  and  the  average  is  about  4'’/i  if  5  =  13  kG.  For  Alcator  C  with  I  = 
5  MeV,  1  =  1.5  MA,  the  local  /3,  is  about  10"/i  while  the  average  fi  is  about  2.5"/'.  In  this  respect  it  is 
useful  to  know  that  conventional  tokamak  plasma  equilibria,  which  are  stable  to  low  m  number  balloon¬ 
ing  modes,  can  be  found  [31]  and  which  have  an  average  p,  of  12‘!<i.  Thus  it  seems  that  the  beam 
plasma  systems  which  we  have  considered  may  well  be  stable  to  ballooning  modes.  Also  stable 
tokamak  plasmas  with  average  )3  of  2.5%  have  been  produced  experimentally. 

Let  us  now  discuss  briefly  a  relevant  experiment  [32].  Intense  electron  beams  have  been  pro¬ 
duced  in  a  controlled  manner  in  Ormak  by  lowering  the  density  to  produce  runaway  discharges.  The 
energy  of  the  beam  electron  is  several  MeV,  the  total  current  is  about  10'^  A,  and  the  confined  beam 
energy  is  more  than  2  kJ.  These  beams  are  stably  confined  in  Ormak  for  35  msec  or  longer.  The 
actual  decay  time  is  about  100  msec,  but  there  is  a  series  of  steps  in  a  current  decrease,  which  shortens 
the  time  to  35  msec.  These  steps  presumably  result  from  some  instability.  Thus  in  Ormak,  an  intense 
relativistic  electron  beam  can  be  confined  for  times  relevant  to  the  scheme  considered  here. 


Vn.  CALCULATIONS  OF  THE  DYNAMICS  OF  TOKAMAK  PLASMAS 

In  this  section  we  solve  Eqs.  (1),  (2),  (3),  (12),  (17),  and  (18)  numerically  for  tokamaks  ranging 
from  small  research  tokamaks  to  large  reactors.  In  all  cases  we  take  r„B/RB.  -  q  We  first  con¬ 
sider  Versator,  a  small  research  tokamak  with  r„  =  7.5  cm,  a  =  15  cm,  /?  =  45  cm,  and  5  =  10  kG. 
A  beam  with  voltage  and  current  of  3  x  10-'  V  and  6  x  10^  A  is  shot  in,  and  the  total  plasma  current  is 
taken  as  equal  to  the  60  kA  beam  current.  Assuming  the  pulse  duration  is  the  time  for  a  beam  ion  to 
go  around  the  torus,  the  beam  time  and  energy  are  t  =  3  x  10"’  sec  and  £  =  6  x  10'  J.  At  f  =  0,  the 
electron  temperature  is  100  eV  and  the  ion  temperature  is  50  eV.  Figure  4(a)  shows  the  time  depen¬ 
dence  of  Tf,  7],  /,  and  =  0)  (/  is  the  total  current)  for  a  proton  beam  shot  into  a  hydrogen 

plasma  with  a  density  2  x  lo".  Even  with  such  a  modest  energy  beam,  there  is  a  large  amount  of  elec¬ 
tron  and  ion  heating.  Figure  4(b)  shows  the  dependence  of  maximum  electron  temperature,  maximum 
ion  temperature,  and  lifetime  on  initial  plasma  density.  The  advantage  of  higher  density  is  apparent. 

Next  we  have  performed  a  series  of  calculations  for  PLT,  with  r„  =  20  cm,  o  =  40  cm,  and  R  = 
140  cm.  We  consider  a  proton  beam  with  V  =  2.1  MeV,  /  =  500  A  (440  kJ)  injected  into  a  hydrogen 
plasma  in  a  magnetic  field  of  40  kG.  The  initial  beam  current  is  the  total  current  and  £,,(0)  =  1  keV 
and  7’(0)  =  500  eV.  In  Fig.  5(a)  is  shown  the  time  dependence  of  total  current,  electron  and  ion  tem¬ 
perature,  and  relative  beam  velocity  as  a  function  of  time  for  plasma  densities  of  10''*.  The  current 
decay  time  is  much  longer  than  the  beam  decay  time  since  the  plasma  rapidly  heats  up,  decreasing  the 
resistivity  so  that  the  current  remains  frozen  in  for  a  long  time.  Note  that  even  though  the  beam 
transmits  forward  momentum  to  the  electrons  they  end  up  going  backward.  The  reason  is  the  induc¬ 
tance  of  the  system.  When  the  electrons  try  to  accelerate  forward,  an  inductive  electric  field  builds  up 
and  drives  them  back.  Figure  5(b)  shows  the  maximum  electron  and  ion  temperature,  discharge  life¬ 
time,  and  the  Q 's  for  a  tritium  beam  shot  into  a  DT  plasma  as  a  function  of  density. 

Here  Q  is  defined  as  the  fusion  energy,  from  both  beam  plasma  and  plasma-plasma  reactions, 
divided  by  the  initial  beam  energy  plus  the  initial  poloidal  field  energy  plus  the  initial  plasma  energy. 
The  a  particle  energy  is  assumed  to  be  deposited  directly  in  the  electrons,  and  the  reaction  rates  are 
those  given  in  the  NRL  plasma  formulary.  Notice  that  temperatures  maximize  at  a  particular  density. 
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This  results  from  the  competition  between  the  electron  energy  confinement  time  increasing  with  den¬ 
sity,  and  the  larger  number  of  particles  sharing  the  beam  energy. 

In  the  remainder  of  this  section,  we  examine  only  thermonuclear  plasmas.  Our  first  study  is  for 
Alcator  C,  with  /■  =  10  cm,  a  =  20  cm,  and  R  =  60  cm.  We  examine  whether  Q~  1  can  be  achieved 
without  pushing  either  the  beam  or  the  tokamak  particularly  hard.  We  examine  plasma  densities 
between  10''’  and  10'^  and  magnetic  fields  between  60  kG  <  <  90  kG.  These  seem  to  be  rather  rou¬ 

tine  operating  parameters  for  Alcator  C.  The  total  current  is  chosen  to  ensure  that  ^  =  1,  that  is  /  =  6 
X  10^  (B/75  kG)  A  and  the  voltage  is  2.1  MeV.  The  reactor  0is  greatly  enhanced  when  we  can  use  a 
beam  current  larger  than  the  total  current,  so  that  there  is  some  current  cancellation  initially. 

First,  we  have  found  that  Q  is  maximized  for  a  tritium  beam  shot  into  a  plasma  which  is  one 
quarter  tritium  and  three  quarters  deuterium.  Figure  6  is  a  plot  of  Q  as  a  function  of  density  for  three 
different  amounts  of  current  cancellation.  In  the  bottom  curve,  the  initial  beam  current  is  the  total 
current;  in  the  middle  curve,  the  initial  beam  current  is  1.5  times  the  total  current;  and  in  the  top  curve 
it  is  double.  For  the  top  curve,  where  B  =  IS  kG,  the  initial  beam  energy  is  780  kJ.  Here  Q  depends 
very  weakly  on  B  but  depends  strongly  on  plasma  density  and  net  beam  current. 


Q 


Fig  6  —  Q  isa  Tunciion  of  densiiy  and  beam  current 
10  lotal  curreni  for  .Mcaior  C 


A  plasma  having  0—1,  while  not  a  pure  fusion  reactor,  is  of  interest  as  a  fission  fusion  hybrid 
[331.  It  may  be  that  something  like  Alcator  C  could  be  used  to  study  this  process  on  a  small  scale.  It  is 
particularly  interesting  that  in  this  case,  neither  the  beam  nor  the  tokamak  is  being  pushed  particularly 
hard,  so  there  is  room  to  scale  up  both  if  for  some  reason  our  calculations  are  overly  optimistic. 
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We  next  examined  whether  a  thermonuclear  ion  temperature  can  be  achieved.  For  Alcator  C,  we 
find  that  an  ion  temperature  of  10  keV  can  be  achieved  with  no  current  cancellation  if  n  =  10'\  Z#  = 
160  kG,  and  the  beam  parameters  are  F  =  5  MeV,  /=  1.5  MA  (beam  energy  of  1.5  MJ).  The 
required  field  can  be  reduced  by  allowing  for  current  cancellation.  For  instance  we  find  that  if  n  = 
10'\  B  =  100  kG,  /  =  1  MA,  /  (Beam)  =  2  MA,  and  F  =  2.1  MeV,  an  ion  temperature  above  10 
keV  can  be  obtained.  For  TFTR  with  r  =  50  cm,  a  =  100  cm,  R  =  250  cm,  «  =  2  x  10’^  or  4  x  lO'"*, 
B  =  50  kG,  and  /  =  2.5  MA,  a  5  MeV  beam  having  7.5  MA  can  reach  the  ignition  temperature.  The 
main  thing  which  allows  for  ignition  is  the  enormous  beam  power.  In  each  case  the  beam  energy  is 
deposited  in  the  plasma  in  about  100  msec,  implying  a  power  of  tens  of  megawatts  for  Alcator  C  and 
hundreds  of  megawatts  for  TFTR. 

VIII.  SUGGESTIONS  FOR  SMALL-SCALE  EXPERIMENTS 

Although  the  injection  schemes  discussed  in  Sections  V  to  Vll  may  be  somewhat  uncertain,  it  is 
very  useful  to  know  that  they  can  be  tested  on  small-scale  experiments  in  either  linear  or  toroidal 
geometry.  Since  cross  field  propagation  has  been  [5,6],  it  remains  to  show  that  the  beam  could  be 
trapped  and  guided  by  a  target  plasma.  Thus  the  beam  could  be  injected  nearly  parallel  to  the  magnetic 
field  into  a  cylindrical  chamber  which  has  a  plasma  localized  in  the  center  of  it  as  shown  in  Fig.  7. 
According  to  the  theory  presented  here,  when  the  beam  collides  with  this  plasma,  the  polarization  field 
should  short  circuit  and  the  beam’s  cross  field  motion  should  stop.  That  is,  the  plasma  should  form  a 
channel  for  the  beam.  A  very  interesting  experiment  then  would  be  to  measure,  and  attempt  to  optim¬ 
ize,  the  efficiency  of  beam  propagation  from  the  diode  to  a  detector  far  down  the  second  cylindrical 
chamber. 


DIODE 


Fig.  7  —  Schematic  of  an  experiment  to  test  injection  into  a  full-density, 
low-volume  plasma  in  linear  geometry 


An  analogous  experiment  should  also  be  possible  on  a  small  tokamak.  The  idea  here  would  be  to 
use  as  a  limiter  a  pin  inserted  from  the  inside.  This  pin  would  define  the  largest  radius  flux  surface 
with  plasma.  The  beam  would  shot  in  from  the  outside.  The  diode  could  be  either  be  in  the  vacuum 
region  of  the  tokamak  or  in  a  guide  tube  outside.  Once  the  beam  hits  the  plasma  it  would  be  trapped 
by  it  as  long  as  Eq.  23  is  satisfied.  Since  the  beam  orbits  are  displaced  outward  from  the  flux  surfaces, 
the  beam  orbits  should  miss  the  pin  limiter  on  the  inside  and  the  beam  plasma  system  should  be  in 
equilibrium.  This  is  shown  schematically  in  Fig.  8. 

It  is  also  possible  to  do  experiments  on  current  generation  in  linear  geometry.  Imagine  a  diode  in 
an  evacuated  guide  tube.  The  diode  produces  a  charge-neutralized,  intense  pulsed  ion  beam  which  pro¬ 
pagates  down  the  guide  tube.  A  short  way  down,  the  beam  comes  to  a  foil  which  separates  the  vacuum 
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from  a  low-density  plasma.  The  decay  of  the  return  current  in  this  plasma  will  measure  how  much 
current  can  be  generated  for  various  plasma  densities.  Such  experiments  would  also  demonstrate  how 
long  it  takes  to  generate  the  current,  what  current-driven  instabilities  are  excited,  etc.  This  sort  of 
experiment  could  give  important  information  as  to  the  background  plasma  density  and  vertical  field 
needed  to  trap  an  ion  beam  in  a  toroidal  chamber. 


Both  of  the  injection  schemes  which  capitalize  on  the  vertical  drift  of  the  ion  beam  rely  on  the 
fact  that  the  beam  is  charge  neutralized  and  that  motion  of  the  unconfined  background  plasma  will  not 
disrupt  the  beam.  These  concepts  could  also  be  tested  on  a  simple  experiment  ir  linear  geometry. 
Imagine  a  diode  in  an  evacuated  cylindrical  tube  in  vacuum.  A  little  way  dow  )  the  tube,  a  foil 
separates  a  region  of  gas  or  low-temperature  plasma  whose  density  can  be  varied.  In  the  plasma  region, 
the  guide  tube  bends  through  for  instance  a  180“  or  270°  turn.  There  might  or  might  not  be  a  vertical 
field.  A  detector  on  the  other  side  of  the  bend  can  test  whether  the  intense  beam  follows  the  single 
particle  drift  orbit.  Specifically  this  would  determine  the  density  requirement  on  the  background  plasma 
for  efficiently  propagating  the  beam  along  the  curved  field. 


Finally  let  us  note  that  injection  of  an  ion  beam  into  a  tokamak  gives  an  ideal  way  to  study  the  /3 
limitation  of  a  tokamak.  According  to  Lovelace  I25],  the  energy  density  of  the  ion  beam  is  interchange¬ 
able  with  plasma  pressure  in  MHD  stability  theory.  Since  the  injection  time  of  the  ion  beam  is  small 
compared  to  any  MHD  instability,  the  plasma  "pressure"  and  magnetic  pressure  can  be  continuously 
varied  by  changing  the  beam  energy  and  toroidal  field.  Thus  one  can  experimentally  study  a  type  of 
high-beta  tokamak,  and  particularly,  can  examine  the  conditions  for  the  onset  of  MHD  ballooning 
instabilities. 


IX.  CONCLUSIONS 


We  find  that  tokamak  heating  and  current  maintenance  with  intense  pulsed  ion  beams  is  an 
extremely  promising  area  for  future  experimental  and  theoretical  studies.  Our  calculations  indicate  that 
one  shot  from  an  ion  beam  at  the  initiation  of  the  discharge  can  be  sufficient  to  reach  ignition  tempera¬ 
ture  in  either  small-volume  or  large-volume  tokamaks.  Since  ion  beams  are  inherentlj  very  efficient, 
and  they  allow  the  elimination  of  the  entire  steady  state  ohmic  heating  circuit,  they  could  most  likely 
give  rise  to  very  economical  reactor  designs.  Finally  there  are  several  interesting  small-scale  physics 
experiments  that  can  be  performed  to  test  the  concepts  developed  in  this  report. 
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